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Abstract 

We present a Scan Method analysis of the allowed region of the 'p -fj plane using the 
latest input measurements of CKM matrix elements, sin 2/3, ^ mixing, ex, o and 7 . 

In this approach, we make no assumptions as to the distribution of theory uncertain¬ 
ties; rather, we scan over the range of plausible theoretical uncertainties and determine 
confidence level contours in the p -p plane. We determine a from branching fraction 
and CP asymmetry measurements of B decays to all light pseudoscalar-pseudoscalar, 
pseudoscalar-vector, vector-vector, and ai-pseudoscalar mesons and determine 7 from 
TT and Dp modes, thereby including correlations between the angles of 
the unitarity triangle. We parameterize the individual decay amplitudes in terms of 
color-allowed tree, color-suppressed tree, gluing penguin, singlet penguin, electroweak 
penguin, as well as W-exchange and W-annihilation amplitudes. Our procedure ac¬ 
counts for all correlations among the fitted CKM parameters (p, p, A and A). The data 
are consistent with the Standard Model with no need for new physics contributions. 
We also examine example wall plots, i.e., projections of sensitive parameters showing 
correlations among them and regions of preferred theoretical parameters. 


1 Introduction 

The phase of the CKM matrix is responsible for CP violation in the Standard Model 
(SM) p. Unitarity relations of the CKM matrix provide an excellent laboratory to test 
this prediction. The relation Vff^Vud + ^ch^cd + ^th^td = 0 is particularly useful since, in the 
Wolfenstein parametrization [2], it represents a triangle in the p — rj plane. For this test, 
many measurements in the B and K systems can be combined. We perform the test with 
the SCAN Method [3], a frequentist technique in which we make no assumptions as to the 
distribution of theory uncertainties of lattice parameters and \ Vuh\ and \ Vch\] rather, we scan 
over them using grid or MG methods. We determine confidence level (CL) contours in the 
p — rj plane and test the hypothesis that the SM is correct using a test compared to the 
alternative that the SM is incorrect. 
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2 Fit Methodology 


We perform baseline fits in which we combine measurements of CKM matrix elements, BB 
mixing, CP violation in the kaon system and angles of the unitarity triangle in the 
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Here, M-p terms represent other measurements such as charm and top quark masses and 
Tn terms represent lattice parameters that have theory uncertainties. We also perform “full 
fits” in which we use as inputs the measurements that determine a and 7 rather than the 
derived values of these angles. Confidence regions are determined by comparing values 
with a critical value instead of looking for a change in The algorithm for a (1 — <ac)[3 
confidence region in d dimensions of a p dimensional parameter space with n measurements 
is as follows: We first determine the acceptance region, at the ac significance level, by 
determining the critical value Xc such that P{x^ > ^ + d\Ho) > cTc, where Hq is the 

hypothesis of the SM, and n — p -\- d is the number of degrees of freedom. The confidence 
region is then given by all those points in the d dimensional parameter subspace for which 
<xl, under Hq. 


3 Fit Results in the p — f] Plane 

In the baseline fits, we fit 23 measurements (|I4d|, |I4s|, \Vcb\, \Vub\, | Vsl, | Vdl, \Vtb\, ^k, ^rna, 
ArUs, sin2/3, B{B riy)jBsJBs/fBd^BBs^BBjBBa,BK,mt,mc,TBa,TB,) listed in 

Tablesand [2[to 13 parameters (p, p. A, A, /s,, mt^rric, tb^^tb,) 

The PDG [4] uses scaling factors of 2.6 and 2.0 for averaging \Vub\ and \Vcb\ results from 
inclusive and exclusive modes, respectively. This procedure increases all errors. Figure [T] 
(left) shows the overlay of 1 — etc confidence level (CL) contours of all accepted fits. 

We also perform “full fits” in which we use 256 branching fraction and CP asymme¬ 
try measurements instead of the a and 7 direct inputs. These fits, with 114 parameters, 
uniquely among the procedures in common use [El US] account for possible correlations 
among ct,/3 and 7 . We currently include all branching fraction and CP asymmetry mea¬ 
surements of B to pseudoscalar-pseudoscalar (PP)^ pseudoscalar-vector (PH), vector-vector 
(VV) and ai-pseudoscalar (aiP) modes to determine a. We include tree, color-suppressed 
tree, penguin, singlet penguin, VF-exchange, VF-annihilation, and EW penguin amplitudes 

*ac is a value around 5% 

Wrom a fit to branching fractions and CP asymmetries of B ^ PP, PV, VV, aiP decays (see below). 

^From a fit to measurements of P ^ decays (see below). 
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Table 1: Measurement inputs. 


Input 

Value 

Ref 

Input 

Value 

Ref 

iKoi 

(4.09 ±0.069 ±0.09) x 10“^ 

in 

mt 

(173.07 ± 0 . 88 )GeV/c^ 

a 

\Vub\ 

(4.15 ±0.31 ±0.39) X 10“3 

B! 

rric 

(1.275 ±0.025)GeV/c2 

m 

IK.I 

0.2252 ±0.0009 

M 

Arrid 

(0.510 ±0.003)ps-i 

m 

\Vud\ 

0.97425 ± 0.00022 

0 

Arris 

(17.761 ±0.022)ps-i 

a 

IKdl 

0.23 ±0.11 

in 


(2.228 ±0.0011) X 10-3 

a 

IKdl 

1.006 ± 0.023 

in 

sin 2/3 

0.682 ± 0.019 

a 

\Vtb\ 

0.97 ±0.08 

in 

a 

85.i7;2° 

t 

B{B Tiy) 

(1.15 ±0.23) X 10“^ 

M 

7 

78.9l?o®3 
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Table 2: Lattice parameters and QCD parameters. 


Lattice 

parameter 

Value 

Ref 

QGD 

parameter 

Value 

Ref 

Jb, 

(228.66 ± 2.0 ± 5.5)MeV 

m 

Vcc 

1.39 ±0.35 

BM 

fBs/fBd 

1.205 ± 0.0086 ± 0.0188 

m 

Vtc 

0.47 ± 0.04 

aim 

Bb. 

1.311 ±0.046 ±0.076 

m 

m 

0.5765 ± 0.0065 

aim 

BbJBb^ 

1.053 ±0.040 ±0.064 

m 

Vb 

0.551 ±0.007 

im 

Bk 

0.7584 ±0.0020 ±0.019 

m 





(up to beyond leading order), as well as SU{3) corrections in defining the amplitudes. 
We parametrize the observables using the Gronau-Rosner approach 

The amplitude ratios and CP asymmetry measurements of and 

B modes determine 7 , while those of B^ B^ T)* 7 r^ and B Dp^ 

modes determine sin(2/3 + 7 ). Figure (right) shows an overlay of 1 — ac CL contours of 
all successful full fits in p — 7 plane. 

Table 3: 1 — ac CL ranges for Unitarity Triangle parameters from baseline fits and full fits. 


Parameter 

Baseline Fit 

without B ^ TV 

Baseline Fit 

with B ^ TV 

Full Fit 

without B ^ TV 

p 

0.069-0.144 

0.073-0.145 

0.081-0.126 

V 

0.320 - 0.395 

0.324 - 0.396 

0.331 - 0.375 

pn 

19.4-24.2 

19.7 - 24.3 

20.1 - 22.8 

an 

79.8 - 90.2 

80.0-90.1 

81.8-88.2 

7[°] 

68.0 - 78.7 

68.0 - 78.2 

70.3 - 77.0 


3 































n 



Figure 1: Contours at 1 — ac CL in the p — rj plane for baseline fits (top) and full fits 
(bottom). 


4 Extraction of a and 7 

We use 181 branching fraction and CP asymmetry measurements of B ^ PP, B PV^ 
B VV, and B aiP modes [6jto fit 94 parameters and determine the a — (3 con¬ 
tour @(1 — ac) CL shown in Fig. (left). The central value of [3 is consistent with the 
measured world average of /3 = (21.51 qy)° extracted from sin2/3 measured in 6 ^ ccs 
modes. Further, we use 56 branching fraction and CP asymmetry measurements of B^ 
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L)(*) 7 r, and D ~modes [ 6 ] analyzed with the GLW [MlEl], ADS [221IM] 
and GGSZ [2l| methods to fit 19 parameters and plot 7-/3 contours @(1 — ac) CL, which 
are shown in Fig. (right). Since the contours depend on we explicitly scan over 

the ratio. The central values for (3 again agree with the world average measured in 6 ^ ccs 
decay modes. 



Figure 2 : The 1 — ac CL contour in the a — f3 plane from fits to branching fractions and CP 
asymmetries of B ^ PP, PV^ VV^ aiP modes (left) and 1 — o^c CL contours in the 7-/3 
plane from fits to p(*) 7 r, and D ~modes (right). The individual 

^ — (3 contours result from a scan over iKi^/V^cfel- The grey-shaded band shows the world 
average of /3 measured in 6 ^ ccs modes. 


Table 4: Determination of the angle a from fits to branching fractions and CP asymmetries 
of P ^ PP, Py, VV, and aiP modes and determination of 7 from rate and CP asymmetry 
measurements of B^ B^ DK*^, B and B^ D~p^ modes. 


Mode 

«[°] 


7[°] 

B ^ PP + PV + VV + aiP 

85.1>? 

21.87;? 


B D^*')h 


22.87;^ 

78.97o®3 


5 Wall Plots 

We construct wall plots to display the correlations among sets of three out of the seven 
parameters with large theory uncertainties. To study the impact of the remaining param¬ 
eters on two parameters displayed, we impose different constraints on the other displayed 
and undisplayed parameters. These studies clearly show that certain regions of the theory 
parameters are favored by the data while others are not. As an example. Figure shows 
correlations for Bb^ versus Bk versus (fop) and Bb^ versus \Vub\ versus /b^ (bottom) 
for 68 % CL (left) and 1 — o^c CL (right). Orthogonal solid lines in each plane show the ±16 
theory error range. Outer black contours result from a probability requirement of >32% or 
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> 5%. Inner black contours result from a ±15 requirement on all undisplayed parameters 
while colored solid contours result from a ±15 requirement on the out-of-plane variable. 
Constraining the out-of-plane variable to its central value yields the colored dashed con¬ 
tours while the black dashed contours result from constraining undisplayed variables to 
their central values. The plots show that larger values of Bk and fss and smaller values of 
\Vub\ are favored. 




Figure 3: Wall plots for Bb^ versus Bk versus (fop) and Bb^ versus \Vub\ versus /b^ 
(bottom) for 68% contours (left) and 1 — ac contours (right). 
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6 Conclusion 


Baseline fits and full fits performed with the SCAN method yield results that are in good 
agreement with the SM without making any assumptions about the distribution of theory 
uncertainties as is done for UT^t [12] and CKMfitter [13] analyses. The full fits how a 
slightly larger allowed region in the 'p — fj plane. They account for correlations among the 
observables cr, [3 and 7 . These correlations, however, are ignored in UTfit and CKMfitter 
analyses. Wall plots allow the determination of whether any potential SM discrepancy 
originates from the values of theory parameters or from experimental measurements. The 
measurements prefer larger values of Bk and and smaller values of \Vub\‘ 
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